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ABSTRACT 

We report the discovery of a substellar companion to 2MASS J02192210-3925225, a young 
M6 7 candidate member of the Tucana-Horologium association (30 — 40 Myr). This L4 7 com¬ 
panion has been discovered with seeing-limited direct imaging observations; at a 4" separation 
(160 AU) and a modest contrast ratio, it joins the very short list of young low-mass companions 
amenable to study without the aid of adaptive optics, enabling its characterization with a much 
wider suite of instruments than is possible for companions uncovered by high-contrast imaging 
surveys. With a mo del-dependent mass of 12-15 Mj^p, it straddles the boundary between the 
planet and brown dwarf mass regimes. We present near-infrared spectroscopy of this companion 
and compare it to various similar objects uncovered in the last few years. The J0219-3925 system 
falls in a sparsely populated part of the host mass versus mass ratio diagram for binaries; the 
dearth of known similar companions may be due to observational biases in previous low-mass 
companion searches. 

Subject headings: 


1. Introduction 


The spectacular discoveries brought on by high- 
contrast imaging of exoplanets around nearb:v 


nearby 

i][2ooi 


stars (e.g., HR8799 ’s system. iMarois et al.l _ 

20101 : B Pictoris b, Lagrange et al. JoiOt Foma- 
Ihaut’s companion, Kalas et al.l 2008 ) eclipse a 
fact that is often overlooked: most of the region 
around a star where a planet would be gravita¬ 
tionally bound for Gyrs is readily accessible to 
seeing-limited observations for nearby stars. A 
handful of planetary-mass objects have indeed 
been uncovered in wide-held surveys, such as GU 
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Psc b and Ross 458(AB)c, respectively 2000 and 
1200 AU from their host s llGoldman et al. l 2010t 
Burningham et all I2OIII : iNaud et al.l 2014 ). and 
the presence of companions on such wide orbits 
presents a signihcant challenge to planetary for¬ 
mation models. 

These companions are on orbits too large 
for in situ formation by either core accretion 
( Pollack et al. 1996 : Alibert et al. 2005 ) or grav- 
i tational instability within protoplane t ary d isks 
( Cameron 1978 : Boss 1997 ). Vorobvov ( 2013 ) es¬ 
tablished that wide companions may form in wide 
orbits within gravitationally unstable protoplan¬ 
etary disks, but only around massive (> 0.7 Mq) 
hosts with massive (> 0.2 Mq) protoplanetary 
disks, and cannot account for the discovery of 
such companions around K stars or later-type 
hosts . Turbulent fragmentatio n of a pre-stellar 
core ( Padoan fc Nordlund 2002h is a viable alter¬ 
native to explain the existence of such systems. 
Outward migration through planet-planet inter¬ 
action is also plausible, but remains to be tested 
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as a credible mechanism at su ch separations and 
mass ratios ( Veras et al. 2009ll . Statistical anal¬ 


yses of the mass, age and separation of planets 
uncovered through high-contrast imaging cam¬ 
paigns suggest that these companions constitute 
a low-mass tail of the brown dwarf distribution 
and are therefore the results of disk or cloud frag¬ 
mentation. The dearth of low-mass (< 5Mj) 
companions found by direct imaging surveys with 
adaptive optics (AO) can be explained in that 
context ( Brandt et al.l[2014 1. This analysis mostly 
contrains the occurence of companions out to 100 - 
200 AU; whether this paucity of companions holds 
for more distant brown dwarfs and planetary-mass 
companions remains to be determined in a statis¬ 
tical framework. Overall, much work needs to be 
done to understand the origin and demographics 
of these distant companions, both theoretically 
and observationally. 

These distant companions have distinct ad¬ 
vantages compared to either isolated planetary- 
mass objects or planets detected through AO sur¬ 
veys. The presence of a host star allows a cross¬ 
calibration of the properties of the planet (paral¬ 
lax, mass, age, metallicity, membership to a young 
group), and the projected separation to that host 
is sufficient to allow direct study through means 
not compatible with extreme-AO, such as accurate 
spectro-photometry, high-resolution spectroscopy, 
optical imaging, etc. Indeed, the detailed analysis 
performed with intermediate-resolution spectro¬ 
photometry on GU Psc b and Ross 458(AB)c will 
be impossible to obtain in the near future for ex¬ 
oplanets uncovered by AO surveys. 

As distant companions of nearby young stars 
provide important benchmarks to understand self- 
luminous gas giants, we undertook various seeing- 
limited observations usii^ the SIMON near- 


infrared spectro-imager ( Dovon et al. 200 n[l at 


the CTIO-I.5m telescope and GMOS-S at Gem¬ 
ini South in order to identify such new objects 
through their distinctive far-red and near-infrared 
colors. We report here a first discovery from 
this survey; a co-moving companion to 2MASS 
J021922I0-3925225 (J02I9-3925), an M 67 candi¬ 
date member of the Tucana-Horologium associa¬ 
tion (THA) that has a low-gravity L4j compan¬ 
ion. 

The survey that led to this discovery is de¬ 
scribed in section [2] The discovery and photo¬ 


metric follow-ups of the companion and its host 
are described in section [H Results are described 
in section |4] and discussed in section [S] 

2. The survey 

J02I9-3925 has been observed as part of a sur¬ 
vey of 300 stars conducted at the GTIO-1.5m 
telescope with the SIMON. This survey was un¬ 
dertaken following the discovery of the planetary- 
mass companion around the M3 AB Doradus 
member GU Psc ( Naud et al 11 ^ . in an attempt 
to identify additional comparable planetary-mass 
companions and assess their overall frequency. 
The target sample consists of both confirmed and 
strong members of young moving groups within 
70 pc (< 120 Myr; (3 Pictoris, AB Doradus, THA, 
Columba, Garina and Argus), and it wa s as- 


sembled from obje cts in Malo et al. ( 20131) and 


iGagne et al. ( 20HJ) . For each star, the Bayesian 
Analysis for Nearby Young Asso ciatioNs II tool 
(BANYAN II; Gagne et al.ll20149) provided a sta¬ 
tistical distance estimate, allowing us to derive a 
projected distance for putative companions. 

By combining S IMON J-band ima ging with 
WISE photometry ( Wright et al. 2010l) . we iden¬ 
tified objects that had a similar position to brown 
dwarfs or known very-low gravity L dwarfs in an 
Mwi vs J — W2 diagram, but that did not have 
optical digitized sky survey (DSS) counterparts. 

In a way similar to the survey that allowed 
the discovery of GU Pscb, candidates were then 
followed-up with Gemini-South using deep i and z 
band imaging to confirm the very red i — z color ex¬ 
pected for an L or T dwarf. The selection methods 
will be detailed in a future paper (Artigau et al. 
in preparation), as refinement of selection methods 
and follow-up of candidates are still ongoing. 

As the SIMON observations have significantly 
better resolution than 2MASS (^I" versus ^2" 
full-width at half maximum) and are significantly 
deeper ( 10 cr at J ~ 18), they provide an oppor¬ 
tunity to identify relatively tight (2 — 6 ") over¬ 
looked companions. We performed a radial profile 
subtraction on all primaries and visually inspected 
the residuals. Through this process, a single can¬ 
didate was identified, at an angular distance of 
about 4" from the M 6 star J0219-3925. No other 
star within our sample presented a similar candi¬ 
date companions within in a 2 — 6 " annulus down 
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to a contrast of A J < 5. 


3. Observation and reduction 


3.1. Imaging 


The J-band discovery dataset of J0219-3925 
was obtained with SIMON on 2013 November 5 
at the CTIO, and follow-up H and itTg-band imag¬ 
ing were obtained on 2014 February 2, with the 
same instrument and telescope. For all 3 imaging 
sequences, we employed a 4-point dither pattern 
along the corners of a 2' x 2' square. At each 
dither position, 15 images (J) and 5 images {H 
and Kg) were taken. All sequences used a 30 s 
per-frame exposure time. The images were sky- 
subtracted using a sky frame constructed from the 
median combination of all science images taken 
on that night within that band. Flat-fielding was 
performed with a flat constructed from images of 
a flat screen. The astrometric solution was per¬ 
formed using a cross-match of the 2MASS cata¬ 
log with field stars. All images were registered 
and median-combined to produce the final science 
frame. 


The contrast ratio between the two components 
was determined by performing a PSF fitting of 
the two J0219-3925 components using two isolated 
bright and nearby field stars (< 2.'5, J > 13.5) as 
input PSFs. The magnitude of J0219-3925B was 
then determined from the 2MASS magnitude of 
its host and its contrast ratio within each near- 
infrared bandpass. The presence of the compan¬ 
ion is unlikely to have significantly affected the 
2MASS magnitudes of its host, as it contributes 
only 2 to 4% of the total integrated flux and is 
resolved 2 FWHM). 


Upon discovery of a possible faint companion to 
J0219-3925, we noticed that images of that field 
usable for science were present in the Gemini sci¬ 
ence archive (See Figure [T]). GMOS-South {i and 
g-band) and Flamingos-II (J-band) spectroscopic 
acquisition images of the host star have been ob¬ 
tained as follow-up obse rvations of the BA NYAN 
All-Sky Survey (BASS; iGagne et al. 2015h . The 
Ai = 6.22 contrast and non-detection in g {Ag > 
5.5) pointed toward an object much redder than 
the mid-M host, prompting for a dedicated spec¬ 
troscopic follow-up of the companion. Inspec¬ 
tion of the A"-band 2MASS images taken in 1999 
shows that the companion is marginally detected 
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Fig. 1.— J0219~3925 system in SIMON, Flamin¬ 
gos 2 (F2), GMOS and 2MASS imagery. Despite 
modest constraints provided by the 2MASS image 
on the position of J0219-3925B, the ^15 yr delay 
between these archival observations and our dis¬ 
covery images provides significant constraints on 
the common proper motion (see section [4.411 . For 
all images. East is left, North is up. The color 
image is a combination of g, i and F2 J-band im¬ 
ages respectively coded as blue, green and red with 
arbitrary scaling. One readily sees that J0219- 
3925 B is much redder than its host. 
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Table 1: Imaging and spectroscopy datasets for the J0219-3925 system 


Instrument 

Date 

Filter 

Number 
of exposures 

Per-frame 
exposure time 

Gomment 
where applicable 

Gemini program ID 

GMOS-S 

17/09/2012 

9 

1 

10.5s 

Spectroscopy acquisition 
GS-2012B-Q-70 

GMOS-S 

02 / 11/2012 

i 

1 

10.5s 

Spectroscopy acquisition 
GS-2012B-Q-70 

SIMON 

05/11/2013 

j 

60 

30 s 

Discovery imaging 

SIMON 

12/02/2014 

H 

20 

30 s 

Photometric follow-up 

SIMON 

12/02/2014 

Ks 

20 

30 s 

Photometric follow-up 

F2 

28/10/2013 

J 

1 

5s 

Spectroscopy acquisition 
GS-2013B-Q-79 

FIRE 

13/02/2014 


2 

306 s 

J0219--3925, spectroscopy 

FIRE 

13/02/2014 


4 

909 s 

J0219-3925B, spectroscopy 

F2 

04/09/2014 

J 

27 

54 s 

Imaging for astrometry 
GS-2014B-Q-72 


and provided valuable constraints on the common 
proper motion of the pair (see Section HlTl and Ta- 
bleE]) . 

3.2. Spectroscopy 


and its companion are compared to field and low- 
gravity objects of comparable spectral types in 
Figure [21 The spectrum of J0219-3925B has an 
average S/N of ~ 30 when sampled at an i? = 800 
resolution. 


On 2014 February 13, we obtained near-infrared 
spectroscopic observations with the Fo l ded-P ort 
Infrared Echelette (FIRE; ISimcoe et al. 2013ll at 
the Magellan 6.5-m telescope. FIRE was used 
in its high-resolution echellette mode, providing 
a spectral resolution R ^ 5000 continuously over 
the entire near-infrared domain (0.8 —2.4^m). We 
used a slit width of 0"6 under a 1"2 seeing and an 
airmass of 1.3 - 1.6 (J0219-3925B) and 1.8 - 1.9 
(J0219-3925A). The reference star for both ob¬ 
jects was observed at the same airmasses. For 
J0219-3925B, we used an ABBA dither pattern 
for improved sky subtraction on this relatively 
faint target, while for the brighter J0219-3925A 
and reference star, the two exposures were taken 
with an AB dither pattern. 


Data reduction was performed using the stan¬ 
dard FIREHOSE pipeline. Flat-held correction 
was performed using hat frames derived from 
dome hat images and telluric absorption at the 
time observation was derived from an AOV star 
(HD 17683) observation taken immediately before 
(J0219-3925B) and after (J0219-3925 A) the sci¬ 
ence integration. The spectrum of the host star 


4. Results 


4.1. Host star properties and membership 

J0219-3925 has been identihed as a member of 
the THA by two teams independently. First, it 
is part of the 129 new late-t ype (K3 and later 
THA members identihed by iKrans et al. 201 
their membership being established by Li and RV 
measurements, but lack parallaxes and is based 
on spectroscopic distances. Spectroscopic htting 
of optical spectrum and overall SED of J0219- 
3925 respectively lead to a spectral type estimate 
of M4.9 ± 1.0 and M5.9 ± 0.3; these values are 
consistent with the near-infrared spectral type of 
M 67 determined here (see section [43]). The radial 
velocity they measure ( 10.6 ± 0.7kms“^) agrees 
within 0.71 kms“^ with the expected radial veloc¬ 
ity of a THA member in that line of sight. This 
is below the internal dispersion of THA velocities 
(^lkms“^) and well below the cutoff in velocity 
difference between members and non-members in 
their analysis (±3kms“^). 

The host star J0219-3925 was identihed in- 
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Table 2: Properties of the J02193925 system 


Short name 

J0219-3925 

J0219-3925B 

2MASS name 

02192210- 

-3925225 

/ia“ (mas/yr) 

107.37 ±2.27 


(mas/yr) 

-34.95 ± 1.65 


Statistical distance^ 

39.4±2.6pc 

Position angle 

173.9 ±0.2° 

Separation 

3.96 ±0.02" 

Sky-plane separation 

156± lOAU 

B 

18.50 


Ag^ 

> 5.5 


6.6 

R 

15.23 


Id 

13.42 ±0.03 


T e 

^2Mass 

11.381 ±0.026 


Jmko^ 

11.321 ±0.026 

15.54 ±0.10 

H2Mass 

10.811 ±0.027 


Hmko^ 

10.855 ±0.027 

14.63 ±0.10 

-^2Mass^ 

10.404 ±0.025 


Kmko^ 

10.444 ±0.025 

13.82 ±0.10 

Wia 

10.148 ±0.023 

W2S 

9.901 ±0.020 

WZa 

9.614 ±0.037 

'^rad 

10.6±0.7« 


V sini 

6.5 ±0.48 


SpT 

M 6 7 ^^ 

L4 7 

EW Ha (A) 

-7.02 


EW LigTos (A) 

639.9 


FeH, 

1.068 ± 0.001 

1.040 ±0.002 

FeHj 

1.07 ±0.01 

1.20 ± 0.02 

Ki 1.169/im 

0.17 ±0.08 

0.6 ±0.4 

Ki 1.177/im 

1.25 ±0.06 

0.71 ±0.32 

Ki 1.244/im 

A 

3.0 ±0.5 

Nai 1.138/im 

2.58 ±0.08 

2.4 ±0.5 

i7-Cont 

0.995 ±0.001 

0.963 ±0.001 

VO, 

1.001 ± 0.001 

1.372 ±0.004 


^Girard et al.|[201ll 

^See lGagne et ani2014cl for a discussion on the definition and the accuracy of the statistical distance. 

‘^Only contrasts were derived in g and z, and not magnitudes measurements as the images were not taken under photometric 
conditions. 

feENIS Consortiumll200'5l 

^Skrutskie et aLll2Q06l 

■^The 2MASS to MKO transform has been determined for J0219-3925 from its FIRE spectrum. The companion’s contrast ratio 
has been determined with SIMON that uses MKO filters. 

^Gutri et al.|[20T3 

^Kraus et al. 1 120141 give a M5.9 spectral type from SED fitting and a spectroscopic spectral type of M4.9. 

*The spectrum of J0219-3925B around the KI 1.244/im is strongly affected by bad pixels on the science array. 
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Fig. 2.— {left) Spectra of J0219-3925, compared to a field M 6 dwarf (LP423) and an M 67 (TWA 8 B). 
Both J0219-3925 and TWA 8 B show an i/-band spectrum that is more peaked that that of the higher gravity 
LP423. The inset shows that full resolution spectra in the 1.12 — 1.20/rm domain; the higher gravity LP423 
shows much stronger Nal absorption at 1.138/im and a stronger KI doublet at ^ 1.17/r m. {right) Spectra 
of J0219-3925B, the youn g companion AB Pic torisb (data from Bonnefov et al.l[2014h and the field L3.5 
dwarf GD 165 B (data from McLean et al.l2003h . The iL-band flux of J0219-3925 is noticeably more peaked 
than the two others, pointing toward a very low gravity (compare with middle panels in Figure [5]). {both) 
Flux within each photometric bandpass has been normalized to that of J0219-3925A or J0219-3925B to 
highlight differences in the shape of the SED within each bandpass; normalization intervals are shown as 
straight black lines above spectra. 


dependently as a candidate member of THA 
as part of the BASS survey. 2MASS and All- 
WISE photometry and astrometry were used 
to derive its membership probability to several 
moving groups in the Solar neighborhood. The 
BANYAN II tool compare the input parameters 
to spatial, kinematic and photometric models of 
young moving groups and the field population 
using a naive Bayesian classi fier. We refer the 
reader to iGaene et al.l ( 2014cll for details on this 
analysis. As J0219-3925 had been identified as 
a high probability candidate of THA, it was in¬ 
cluded in our search for distant companions before 
the publication of the bulk of the survey. The 
BANYAN II tool gives a P=99.34% probability 
that J0219-3925 is a member of THA. When such 
measurements are available, the BANYAN H tool 
can use the radial velocity and/or parallax to con¬ 
strain membership probability further. Using the 
radial veloci t y of 10.6 ± 0.7kms“^ measured by 
iKraus et al. ( 2014 b the THA membership proba¬ 
bility of the host star increases to P=99.94%. 

We can set additional constraints on the age of 
this object since it displays lithium (jKraus et al 


20141) and typical signs of low-gravity in its NIR 
spectrum (see left panel of figure [2] and Sec¬ 
tion Within that spectral type, the presence 
of lithium (EW= 639.9 mA) points to an age be¬ 
low ^ 125 Myr. Adding this age constraint in the 
BANYAN H analysis raises the THA membership 
probability further (P=99.93%). In addition to 
a membership likelihood, the BANYAN analysis 
provides a kinematic distance to an average pre¬ 
cision of ^ 10 % for young moving group members 
(See Figure 5 iniMalo et al.ll2013l and Figure 8 in 


Gagne et al.ll2014(J) . The kinematic distance for 
J0219-3925 is 39.4 ± 2.6 pc. The 7% uncertainty 
includes both the contribution from the uncertain¬ 
ties on the proper motion of J0219 and the space 
velocity scatter of THA members. 

As an additional confirmation, we used J0219- 
3925 apparent Ic and J magnitudes, the amplitude 
of the proper motion and the sky position to ver¬ 
ify its m embership probab ility in the BANYAN I 
analysis ( Malo et al]l2013l) . It yields an a priori 
probability of P=96.3% of membership in THA. 
Similarly to the BANYAN H analysis, it is possi¬ 
ble to include the radial velocity as an additional 
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constraint in the analysis. When we do so, the 
probability is increased to Pj,=99.9%. The statis¬ 
tical distance inferred by this analysis is the same 
as predicted by BANYAN II. 

4.2. The Age of THA 

The Tucana and Horologi um Associations were 
disco vere d independently by IZuckerman &: WebbI 
2 OOOI and Torres et al.ll2000l Further studies 


con¬ 


firmed that what was initially considered as two 
associations were more likely to constitute a sin¬ 
gle young moving group to be c alled the Tucana- 
Horo logium association (THA; IZuckerman et ah 
I2OOIII . THA is part of the so-called Great 


Aust ral Young Association (GAYA; [Torres et al 


2008!) . including the Colum ba and Garina as- 
sociations. Severa l stud ies ([Zuckerrnan fc Sond 


2004 : [Torres et al.l 2008 : Kiss et ahT 2011 ) pro¬ 
posed members of THA based on their galac¬ 
tic space velocity, galactic position, and signs 
of youth, yielding a sample of members with 
spectral type between AIV and early-M spawn¬ 
ing a d istance range of 36 to 71 pc. Recent 
studies ([Rodriguez et al. 2013: Moor et al. 20131: 


Malo et al.l201^ Gagne et al ] [2014ai bl[d: [Kraus et al 


2014[) proposed more than 200 strong THA can¬ 
didate members, which still require more robust 
kinematics to confirm their membership. 

The age range of THA was estimated to be 10 — 
40 Myr ( Zuckerrnan et al.[[2001 ; Zuckerrnan fc Webb 
2000 [) based on various age indi cators (Ha, 
lithiu m, HR diagram). More recently, [Kraus et al.[ 
2014[ derived an average isochronal age of 30 Myr 
using a new sample of 142 candi date members 
combined with the BGAH models (Baraffe et al.[ 


19981) . 

In addition to the isochronal age, the lithium 
depletion boundary (LDB) is a key indicator to 
determine the age of the low-mass star popu- 
lation. This method was used by [Kraus et al 


2014[ to determine the THA lithium depletion 
age of 41±2 and 38±2 us ing the BGAH and 


D’Antona fc Mazzitellil 1997 models, respectively. 


This discrepancy between isochronal and LDB 


(Song et al.[[2002[:[Yee & JenseiJ201C 

: Binks & Jeffried 

2014). Recentlv. Malo et al. 2014a 

have shown 


that usi ng new Dartmouth Magne tic evolutionary 


age is in better agreement with the LDB age for 
young low-mass stars. In general, isochronal ages 
are revised upward with t he inclusion of mag - 
netic fields. As shown in [Malo et ^ ()2014b[) . 
magnetic field strength measurements from high- 
resolution spectroscopy would further constrain 
its age. Pending such measurement, we conser¬ 
vatively adopt an age range of 30 — 40 Myr for 
the THA association, and thus, for J0219--3925 
system. 

4.3. Spectral type and Gravity Indicators 

We used the method of K. Cr uz et al. (in prepa¬ 
ration; see Cruz &: Nuhed 2007 ) to assign a spec¬ 
tral type to both components of J0219-3925. This 
method consists of a visual comparison with a grid 
of spectral templates while normalizing each NIR 
band individually. Inspecting the slope and shape 
of several features in each band allows to choose a 
template that best matches the observations. Our 
sequence of field, intermediate-gravity and low- 
gravity templates were build by median-combining 
various objects within each spectral types and 
gravity class. The spectr a used to build the tem- 
plates were obtained from [Allers fc Liu[ ([20131) and 
the SpeX Prism Spectral Librarie^. Both the vi¬ 
sual comparisons of J0219-3015 and b yielded best 
matches to very low gravity templates; J0219- 
3015 and its companion were assigned a spectral 
type of M 6 7 and L4 7 respectively. The gr avity 
classification scheme of Allers fc Liul ()2013l ) was 
subsequently used to confirm that both objects 
have weaker alkali lines compared to field dwarfs 
of the same spectral types, and that J0219-3015B 
has stronger VO absorption and a triangular¬ 
shaped iL-band continuum (which was already ap¬ 
parent from the visual comparison). Both objects 
were consistently categorized as very low gravity 
dwarfs by this index-based classification scheme. 
We show that they both display typical spectro¬ 
scopic signatures of a low-gravity, including lower- 
than-normal alkali (K I, Na I and FeH) equiv¬ 
alent widths. The lower gravity causes a lower 
pressure in the atmosphere, which decreases both 
the effects of pressure broadening (responsible for 
the lower alkali equivalent widths) and collision- 
induced absorption (GIA) of the H 2 molecule. The 
triangular H-band continuum is shaped by water 


models (jFeiden fc Ghaboverl[2ni3l ). the isochronal 


^ http://pono.ucsd.edu/-adam/browndwarf s/spexprism 
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absorption. In the case of field brown dwarfs, the 
H 2 CIA redirects part of the flux to the bluer side, 
which masks t he triangular sha pe of the iJ-band 
continuum (see lRice et al. 20101 for further detail). 
Figure [3] illustrates the g ravity-sensitive spec tro- 
scopic indices defined bv lAllers &: Liul ( 2013ll for 
both components of the J0219-3925 system. 

The intermediate gravity and very low gravity 
classifications generally correspond to the 0 and 7 
gravi t y classifications in the optical ( KnIcDatrick 
2005 : Kirkpatrick et al.l 20061 : Cruz et alT 2009ll . 


hence we choose here to use the Creek-letter 
nomenclature even though our gravity classifica¬ 
tion was done in the NIR. 


4.4. Kinematics and Common Proper Mo¬ 
tion 

Astrometric measurements were performed on 
GMOS-S f-band, F2 J-band and SIMON imaging 
obtained on 2013 November 5. SIMON imaging 
obtained in early 2014 were not included as they 
were taken under poor seeing conditions 2") 
and do not constrain proper motion. Archival 
2MASS imaging were also used; it provides only 
modest constraints on the PA and separation, but 
the long time baseline (^15 years) makes it a use¬ 
ful additions to confirm common proper motion. 
Astrometric errors were estimated to be at the 
^0"03 level for F2 and GMOS-S data and ^O'.'l for 
SIMON imaging. Astrometric measurements are 
shown in Figure |4l The for the co-moving case 
is 5.9 for 8 degrees of freedom. The co-moving case 
is equivalent to a 0.9 —cr event in a Gaussian distri¬ 
bution while the background object model would 
correspond to a 6.0 — a event. While this demon¬ 
strates that J0219-3925B is comoving to within 
astrometric uncertainty, any interlopers would be 
an L dwarf at roughly the distance of J0219-3925 
with a significant proper motion that could, con¬ 
ceivably, match that of J0219-3925 within astro- 
metric uncertainties. The strongest argument in 
favour of the two objects forming a physical pair 
arrises from the relative rarity of field L dwarfs per 
sky surface unit. 

Finding a low-gravity L dwarf within 4" of an 
unrelated young M star is very unlikely. Gonsid- 
ering that we started with a sample of 300 young 
stars, the area within a 2 — 6" annulus around 
these stars covers 5.6 x I0“® of the entire celes¬ 
tial sphere. With an L dwarf spatial density of 


-3.8 X 10-3 pc-3 


(jGruz et al.ll2007ll . there should 
be ^5500 L dwarfs within 70 pc. The likelihood 
of finding a coincident L dwarf to a young star in 
our sample is about ~ 3 x 10“'^ without consider¬ 
ing the fact that the L dwarf itself also displays a 
low surface gravity. 

Orbital motion of the pair is expected to be 
close to the detection threshold. With a total 
mass of ^0.12 Mq (see table [3]) and a separa¬ 
tion of '^lOOAU, the orbital period is on the or¬ 
der of 6000yr. Assuming a face-on orbit, this 
leads to an orbital motion of 4mas/yr. Orbital 
velocity is on the order of 800 m/s and within 
reach of upcoming high-resolution i nfrared spec¬ 


trogr aphs such as VLT/CRIRES-I- (|Follert et al 


20141) . Both measurements are challenging but 


possible with existing facilities, but not with the 
dataset in hand. Proper constraints on the dif¬ 
ferential motion within the system will allow one 
to constrain whether the orbit of J0219-3925B is 
consistent with a circular orbit or a highly ec¬ 
centric one. Gonstraints on orbital eccentricity 
for a few systems similar to J0219--3925 could set 
strong limits on the plausible formation mecha¬ 
nisms. These measurements would provide con¬ 
straints on the J0219-3925 system similar to the 
ones obtained for a few simi lar companions in 
the planetary-mass regime by iGinski et al. 2014i 
These authors hnd that fitted orbits suggest that 
distant companions’ motion is consistent with ec¬ 
centric orbits; whether this holds for the J0219- 
3925 system remains to be determined. 


4.5. Model Fitting with Spectroscopy and 
Photometry 

To estimate fundamental parameters of mass, 
effective temperature, gravity, and radius for both 
J02I9-3925 and J02I9-3925B we used two dis¬ 
tinct approaches with the same set of theoretical 
models. First, we constrained the bulk properties 
of J02I9-3925A and J0219-3925B (mass, effec¬ 
tive temperature, surface gravity, radius and lumi¬ 
nosity) by comparing the absolute magnitudes in 
each near-infrared bandpass with BT-Sett l mod- 
el^l (CIFIST 2011 opacities; lAllardll2014ll at 30 
and 40 Myr. Absolute magnitudes were estimated 
from the kinematic distance; uncertainties in the 
kinematic distance were assumed to be 7% (see 


^https: //phoenix. ens-lyon.fr/Grids/BT-Settl/CIFIST2011/ 
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Fig. 3.— Spectral indices as defined bv lAllers fc Liu 20131 for J0219-3925 and J0219-3925B, intermediate 
(INT-G or /3) and very-low gravity (VL-G or 7 ) M and L dwarfs, the field sequence (thick, blue line) and its 
scatter (blue shaded region). The dotted line represents the delimitation between INT-G and VL-G regimes. 
Spectral types were offset by small (< 0.15) random subtypes so that vertical error bars can be distinguished. 
For all indices, J0219-3925B falls outside the enveloppe of field objects, clearly highlighting its low gravity. 
Four indices track both metallicity and surface gravity trends (KIj, FeHj, iJ-cont and FeH^); J0219-3925B 
does not show clear signs of peculiar metallicity, following the general trends of very-low gravity members. 
The M 67 J0219-3925 falls in a part of the diagram where NIR low-gravity indices are less efficient than for 
its L 47 companion. Alkali line equivalent widths are systematical lv lower than field obje cts and consistent 


with very-low grayity objects. Indiyidual objects were drawn from lAllers fc Liull2ni3l and iManiayacas et al 


2014. 
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Fig. 4.— Position angle and separation for the 
2MASS (1999 July 27), CMOS i-band (2012 
Novembre 2), F2 J-band (2013 Octobre 28, 2014 
Septembre 4) and SIMON J-band (2013 Novem¬ 
bre 5) imaging. PA and separation are constant 
to within astrometric uncertainties. Insets in each 
panel shows the astrometric measurements ob¬ 
tained since 2012. Dashed line shows the mean PA 
and separation, while the continuous line shows 
the expected values for a distant background ob¬ 
ject. 


subsection 14.1|) . leading to a 0.14 mag uncertainty 
on the absolute magnitudes. Table [3] compiles the 
value derived for each photometric bandpass, and 
the mean value for all three bandpasses is used as 
the best estimate. Upper and lower bounds are av¬ 
eraged for this best estimate, but not divided by 
the square root of the number of measurements 
as these are correlated (i.e., same age assumption, 
same uncertainty on distance). Overall, the mass 
estimate range of 12-15 Mj^p sets the companion 
at the brown dwarf/planet regime limit. Uncer¬ 
tainties for log g values are small (0.02 and 0.06 
respectively) as the surface gravity varies very lit¬ 
tle with effective temperature for these relevant 
masses and the main contribution to this uncer¬ 
tainty arrises from the uncertainty on the age of 
the system. The relative uncertainty on the mass 
ratio (q) is smaller (^ 6 %) than the relative un¬ 
certainty on either component’s mass (^ 10%) as 
the derived masses for both components correlate 
through the distance and age estimates. 

The second approach to constrain the proper¬ 
ties of J0219-3925B has been to perform a 
fit between the observed and theoretical spectra 
in 3 individual bandpasses (U -I- J, H and AA; 
respectively 1.00 — 1.33/xm, 1.45 — 1.81/iin and 
1.95 — 2.40^m). For each bandpass, spectra were 
normalized over part of the domain (horizontal 
lines in Figure [S|). By performing this normal¬ 
ization, we assume no prior knowledge of distance 
and absolute magnitude. The spectral fitting anal¬ 
ysis has been performed by steps of 0.5 dex in 
logg between 3.5 and 5.5, and over the 1500 K 
to 1800 K domain. Models at intermediate grav¬ 
ity (4.5) better fit the SED. Within the H band 
(middle columns), the higher-gravity (5.5) models 
show the distinctive flattening of the SED, which 
is the hallmark of old, field, L dwarfs. Within 
the K band (right), there is a clear shift of the 
SED peak from ~ 2.1/xm to ~ 2.25 ^m between 
logg = 5.5 and logg = 3.5, for all temperatures, 
with 2M0219-3925B being intermediate between 
these scenarios. For Y + J and AT, the best-fitting 
model has a temperature of 1700K and logg = 
4.5. The value derived for H is only marginally dif¬ 
ferent with Tes = 1600K and logg = 4.0. These 
values are in very good agreement with those de¬ 
rived from photometry alone (see Table [S]) , with 
Teff = 1683 ± 43 K and log g = 4.24±0.04. 

Overall, the mass of J0219-3925B falls squarely 
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Table 3: Constraints on properties from models. 



J0219-3925A 

J0219-3925B 

Mass (Mj) 

114±13Mjup 

13.0±0.7Mjup 

Mass (M//) 

108±llMjup 

13.9 ± 1.1 Mjup 

Mass (Mic) 

115±13Mjup 

14.8 ± 1.6 Mjup 

Mass (mean) 

113±12Mjup 

13.9±l.lMjup 

q (Mj) 

0.113 ±0.007 

q {Mh) 

0.128 ±0.006 

q (Mx) 

0.127 ±0.007 

q (mean) 

0.123 ±0.006 

Teff (Mj) 

3070 ± 73 K 

1615 ± 41K 

Teff (Mh) 

3047 ± 84 K 

1686 ± 43 K 

Teff (Mk) 

3074 ± 72 K 

1746±49K 

Teff (mean) 

3064 ± 76 K 

1683 ± 43 K 

log Luminosity (Mj) 

-2.22 ± 0.06 Lq 

-3.92 ± 0.03 Lq 

log Luminosity (Mh) 

-2.26 ± 0.06 L0 

-3.84 ± 0.05 Lq 

log Luminosity (M^f) 

-2.22 ± 0.06 Lq 

-3.76 ± 0.06 Lq 

log Luminosity (mean) 

-2.23 ± 0.06 Lq 

-3.84 ± 0.05 Lq 

logg (Mj) 

4.59 ± 0.06 cms“^ 

4.23 ±0.02 cm s-^ 

logg {Mh) 

4.59 ± 0.06cms“^ 

4.24 ± 0.04 cms“^ 

logg (Mx) 

4.59 ± 0.06cms“^ 

4.25 ± 0.06 cms“^ 

logg (mean) 

4.59 ± 0.06 cms“^ 

4.24 ± 0.04cms”^ 

Deuterium (Mj) 

0 

0.77 ±0.06 

Deuterium (Mh) 

0 

0.68 ±0.07 

Deuterium (Mjf) 

0 

0.61 ±0.09 

Deuterium (mean) 

0 

0.69 ±0.07 

Radius (Mj) 

2.75 ± 0.14 Rjup 

1.41 ± 0.02 Rjup 

Radius (Mjj) 

2.68 ± 0.12 Rjup 

1.44 ± 0.04 Rjup 

Radius (M^c) 

2.76 ± 0.14 Rjup 

1.46 ± 0.04 Rjup 

Radius (mean) 

2.73 ± 0.13 Rjup 

1.44 ± 0.03 Rjup 


at the upper limit of the International Astro¬ 
nomical Union definition for an exoplanet - con¬ 
sisting of an object with a true mass below the 
limiting mass for thermonuclear fusion of deu¬ 
terium - and is in orbit around a star, indeed 
J0219-3925 has a mass estimate above the hy¬ 
drogen burning limit for plausible ages of THA. 
Evolution models predict that J0219-3925B will 
have retained ^70 % of its initial deuterium by 
40 Myr. As it is expected to have partially burned 
its deuterium, and considering its brightness, it 
would constitue a good target to spectroscopically 
test the brown dwarf/planet boundary, that is 
expected to be mod erately metallicity-dependent 
( Spiegel et all l201lh and its exact predicted lo¬ 


cation varies by depending on atmo¬ 

sphere models ( Saumon fc Marle^l2008ll . The ex¬ 
act nomenclature for such an object is still to be 
properly defined, and it joins a number of other ob- 
i ects such as J()103-5 5 15(AB~)b and J0122 -2439B 
( Delorme et al. 20131: Bowler et al. 2013 ) at the 
upper-limit of planethood. As the J0219-3925B 
is expected to have burned some of its deuterium 
and formally lies just above the deuterium-burning 
limit, we use the ”B” designation for stellar and 
brown dwarf companions rather than the ”b” in 
usage for planetary companions even though it is 
predicted to be slightly less massive than objects 
initially described as planetary companions (e.g., 
AB Pictorisb). 
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Fig. 5.— Comparison between 2M0219-3925B (black) and BT-Settl models (red, green and blue) in the 
Y+J,H and K spectral regions. To highlight shape differences in SED of individual photometric bandpasses, 
we show all models normalized to a common wavelength interval (shown as a thick horizontal line in each 
plot). For the sake of clarity, only \ogg values of 3.5, 4.5 and 5.5 and temperatures from TeS = 1600K to 
1800 K are plotted. 


5. DISCUSSION 


Figure |6] illustrates the position of J0219- 
3925 B in the J — K color-magnitude diagram 
compared to field ultracool dwarfs and planetary 
mass objects. J0219-3925B follows the overall 
trend of planetary-mass companions in being red¬ 
der by about 0.5 mag than field ultracool dwarfs 
of similar Mfc- 

With a mass ratio q = 0.122 ± 0.006, the 
J0219-3925 system is similar to a heavier version 


of 2MASSW J1207334-393254 (Chauyinetal 


2005af l. consisting of a 4 — 6 Mjup planet around a 
25 —45Mjup brown dwarf (mass ratio 0.13 — 0.16). 
With a separation of > 55 AU, J1207-3932 system 
comes very close to J0219-3925 in terms of bind¬ 


ing energy, with a total mass ^3 times larger 
(^ 40Mjup versus 120Mjup), but a physical 
separation that is also about three times larger 
(>55 versus >160AU). It also comes close to 
the J0103-5515(AB)b system, which consists of 
a 12 — 14M.Tiip com panion to a binary mid-M 
( Delorme et al. l2013ll . We can rule-out the pos¬ 
sibility that J0219-3925 is a triple analogous to 
J0103-5515(AB)b. We only have seeing-limited 
observations of the system, so a hierarchical sys¬ 
tem with an inner component tighter than ^075 
would not be resolved with the dataset at hands. 
The best constraints on the presence of an in¬ 
ner binary comes from the Mj, Mh and 
derived from absolute magnitude versus spectral 
type relations for young field dwarfs. These rela- 
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Fig. 6 .— Color-magnitude plot of field dwarfs, low-gravity and directly imaged companions. J0219-3925B 
is redder than typical field L dwarfs, and falls within the sequence of planetary-mass companions. 


tions predict contrast ratios respectively of 4.36, 
3.78, 3.31 mag (Gagne et ah, submitted to ApJ) 
for single M 67 and L4j components, and these 
values are in close agreement with the observed 
contrasts of AJ = 4.2 ± 0.1, AH = 3.8 ± 0.1 
and AKs = 3.4 ± 0.1. We can therefore conclude 
that neither J0219-3925A nor J0219-3925B is an 
equal luminosity binary. We cannot rule-out that 
either object is itself a high-contrast (> 1 mag) 
binary. 


J0219-3925B als o shares various sim ilarities 
with AB Pictorisb ( Chanvin et al. EoOhbl l. Both 
objects are companions to Tucana-Horologium 
members and therefore share a common age. 
AB Pictorisb as a slightly earlier spectral type 
(L0. 5 ± 0.5) and ^0.4 mag brighter in both J and 
Kg ( Biller et al. 2013ll . The brightness difference 
between the two is consistent with what would be 


expect considerin g their respective s pectral types 
(see Figure 17 in iBiller et al. 2013ll . As shown 


in Figure [21 right panel, J0219-3925B displays 
slightly deeper water bands, especially around 
1.5/un and 1.95/un and a deeper CO bandhead 
longward of 2.29/un, consistent with a later spec¬ 
tral type. This brightness difference corresponds 
to a mass difference of ^IMjup between the two 
objects and a comparative study could be used 
to better constrain the deuterium burning limit 
through high resolution spectroscopy. 

Figure [7] shows the host mass versus mass ratio 
for imaged planetary and low-mass brown dwarf 
companions. Interestingly, the J0219-3925 sys¬ 
tem falls in a relative gap among known systems 
in that diagram. Known systems hosting a low- 
mass secondary (< 40Mjup) either have near¬ 
equal mass; many of these systems have been 
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identified as near-equal luminosity field brown 
dwarfs or have a much smaller q values (typically 
q < 2%) and have been uncovered with other 
observing techniques. Efficient searches for near¬ 
equal luminosity sub-stellar bin aries have been 


large ly performed with HST ('e.g.- lBurgasser et al 


2006l l and laser guide star adaptive-optics imag¬ 


ing (e.g., iLiu et all 2012[) while low-q binaries 
have largely been identified by other techniques 
such as microlensing, transit or radial velocity. 
Wide-held imaging surveys have also uncovered a 
few tens of wide c ompanions to held stars (e.g., 
Baron et al.l l2015l IPeacon et al. 2014 and refer¬ 


ences compiled therein); many of these have mass 
ratios comparable to that of J0219-3925, but 
with signihcantly more massive components. The 
two companions uncovered that come closest to 
J0219-3925 are LP 261-75 and Wolf 940. LP261- 
75 AB consists of a young M4.5 /L6 binary at a 


projected separation of 450 AU (jBurgasser et al 


200,51: iReid fc Walkowicd l2006ll . LP 261-75 does 


not match the kinematic properties of young mov¬ 
ing groups included in the BANYAN II tool, but 
its primary has chromospheric activity levels con¬ 
sistent with an age of 100-200 Myr. This youth of 
the system leads to a low mass estimate for the 
secondary of 15 — 30Mjup. Wolf 940 AB consists of 
an M4-I-T8.5 binary with a projected separation of 
400 AU. The system has an age of 3.5-6 Gyr as de¬ 
rived from the chromospheric activity of the host 
star, leading to a mass estimate of 20-32Mjup, 
making it a close match the LP 261-75 system 
albeit at an older age. 


While pairs similar to J0219-3925 may be in¬ 
herently rare, a relatively straightforward obser¬ 
vation bias could explain the fact that similar 
systems have been overlooked. Hosts with spec¬ 
tral types as late as that of J0219-3925A are 
in general too faint at optical wavelengths to 
close an adaptive-optics loop on, and are there¬ 
fore under-represented in planet-search surveys. 
Furthermore, until very recently, very few nearby 
(< 50 pc) young stars that late were known, and 
most comparable objects had been found in sig¬ 
nificantly more distant open clusters, making the 
search of faint companions even more challeng¬ 
ing. While companions similar to J0219-3925B 
may exist around very nearby M dwarfs, the con¬ 
trast ratio increases significantly with age. At 
an age of I Gyr the pair would have a contrast 


ratio of AJ ^ 8 and an effective temperature 
of ~ 600 K while at 5 Gyr, it will have a con¬ 
trast of A J ^ 10 and an effect ive temperature 
of ^ 370 K ( Beichman et al.l 20I0ll . At these mag¬ 
nitudes and temperatures, J0219-3925B will re- 
spectively have spectral types of T8.5 and YO 
( Dupuv fc LhJ 2012 : Marsh et al. 2013 1. Such 
companions at separations of a few arcseconds 
around nearby mid-M dwarfs could easily have 
been overloo ked by current brow n dwarf searches 
using WISE ( Wright et al. 2010ll or seeing-limited 
near-infrared observations due to the far wings of 
the central star’s point-spread function. gaiaI 
astrometric measurements will not constrain the 
occurence of similar systems in the solar neigh¬ 
bourhood; for a face-on circular orbit and a dis¬ 
tance of 10 pc, the host star in a system simi¬ 
lar to J0219-3925 will display an astrometric ac¬ 
celeration of ^ 1.8//as/yr^, below the detection 
threshold for this mission. The GAIA mission may 
nevertheless uncover similar systems in young as¬ 
sociations through their common proper motion. 
The best prospect to find Gyr-old siblings of the 
J0219-3925 system is with the use of very deep 
near-infrared imaging, either under good seeing 
and proper point-spread function subtraction, or 
with laser-guide star adaptive optics. 

Weak constraints can be set on the occurence 
of companions similar to J0219-3925B around M 
dwarfs for the range of separations explored here 
(100— 5000 AU). The PALMS survey ( Bowler et al. 
2 OI 5 II provides the largest sample to assess the 
abundance of analogs to J0219--3925B, although 
the separation range probed with their adaptive 
optics observations only partially overlaps with 
the range of interest here. They set a < 10% upper 
limit on the occurence of planets between 6 and 
200 AU and < 50% limit for 1.8 - 570 AU (COND 
atmosphere models, circular orbits and 95% con¬ 
fidence level). The discovery of a single object at 
the planetary/brown dwarf limit in our sample is 
therefore well within the constraints set by the 
PALMS survey. Our survey included stars with 
different levels of confirmation regarding their age 
drawn from different sources, so deriving clean 
statistical constraints on the abundance of com¬ 
panions similar to J0219-3925B is non-trivial and 
beyond the scope of this paper. J0219-3925 was 


^http://sci.esa.int/gaia/ 
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drawn from the BASS sample of late-type objects 
(>M4). We observed 67 high-probability BASS 
can didates and 13 low -probability ones (see de¬ 
tails [GagnJ^T^OliQlB); the absence of any further 
candidate in the 2 — 6" separation range suggests 
that objects similar to J0219-3925B are relatively 
rare around M dwarfs, with an occurence rate of 
< 6.8% at the 95% confidence level. 
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